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Abstract: We propose a new Ising spin glass model on Z d of Edwards- Anderson type, 
but with highly disordered coupling magnitudes, in which a greedy algorithm for producing 
ground states is exact. We find that the procedure for determining (infinite volume) ground 
states for this model can be related to invasion percolation with the number of ground 
states identified as 2^, where M = Af(d) is the number of distinct global components in the 
"invasion forest". We prove that M{d) = oo if the invasion connectivity function is square 
summable. We argue that the critical dimension separating M = 1 and M = oo is d c = 8. 
When Af(d) = oo, we consider free or periodic boundary conditions on cubes of side length 
L and show that frustration leads to chaotic L dependence with all pairs of ground states 
occuring as subsequence limits. We briefly discuss applications of our results to random walk 
problems on rugged landscapes. 



1 Introduction 

Theoretical work on spin glasses has focused overwhelmingly on the short-ranged Edwards- 



Anderson (EA) model |10[ and its infinite-ranged counterpart, the Sherrington-Kirkpatrick 



(SK) model P7|j. The EA Hamiltonian is: 



n = - £ J- 



X y0~ x Oy 



<xy> 



where the sum runs only over nearest-neighbor pairs of sites on a regular <i-dimensional 
lattice. We will confine ourselves in this paper to Ising models, i.e., the spins a x take on 
the values ±1. The couplings J xy are i.i.d. random variables chosen from a distribution 
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symmetric about zero. One standard choice is a mean zero Gaussian, though one needn't be 
confined to that case. 

There are a number of open questions [[| associated with the EA model. These include 
the question of whether it displays a thermodynamic phase transition, and if so, in which 
dimensions; the nature of its large-scale dynamical behavior for standard local spin-flip 
dynamics; the relationship between the behavior of this model and others, such as the SK 
model or long-ranged models; and others. The question we are concerned with here is the 
number of ground state pairs the model displays in d dimensions. (Because of the spin- 
flip symmetry in the Hamiltonian, ground states always come in pairs.) To avoid local 
degeneracies, we will assume that as in the Gaussian case, the common distribution of the 
J X yS is continuous. 

Parisi's analysis |23| of the SK model displays the striking feature of an infinite number of 
low-temperature phases [[25] organized in an ultrametric fashion [18|. A number of workers 



in the field have assumed that a similar result applies also to short-ranged spin glasses, and 
to other models with frustration and quenched disorder. 

However, an alternative point of view arose in the mid-1980's following a scaling ansatz 
due to MacMillan |L7| , Bray and Moore || , and Fisher and Huse [12| . Their work led to the 



opposite conclusion ||14j| , namely, that the EA model has only a single pair of ground states 
(or of low-temperature pure states) in all finite dimensions. Thus, over the past decade 
two opposing viewpoints have dominated the literature. The first, following the work of the 
Parisi school, asserts that short-ranged spin glasses have an infinite number of ground states 
(in infinite volume) for all "nontrivial" dimensions, i.e., d > 2. The second, proceeding from 
scaling arguments, argues that such spin glass models possess only a single pair of ground 
states in all finite dimensions. (Both of the above statements are normally interpreted to 
hold for almost every coupling realization.) 

In a recent paper pD[ , hereafter referred to as I, we invented a model (see also ||) simpler 
than (but related to) the EA spin glass in which we could explore these issues more fully. 
The model is analytically tractable, yet complex enough so that it displays a surprisingly 
rich ground state behavior. We will argue that our model undergoes a transition in ground 
state multiplicity at eight dimensions: below it has a single pair of ground states, while above 
it has infinitely (in fact, uncountably) many. 

We are able to demonstrate explicitly the mechanism by which multiple ground state 
pairs arise. This is the first such demonstration of which we are aware for a short-ranged 
spin glass model in finite dimension. The nature of the mechanism is related to a mapping 
of the ground state structure of our model to invasion percolation |IB, [7], p9| . We will show 
that solving the problem of ground state multiplicity in our model requires the solution of 
an interesting problem in invasion percolation. 

To summarize, we propose and analyze a tractable short-ranged finite dimensional spin 
glass model (although, like the SK model, it is not physically realistic). Our techniques 
enable us to study the relationship within the model among quenched disorder, frustration, 
dimensionality, and ground state multiplicity. We will find in particular that the role of 
frustration is subtle but important, and our model is instructive in illustrating the larger 
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role frustration may play in more realistic spin glass models. 

The basic features of our model and its analysis have appeared in I. In this paper, we 
supply a number of arguments omitted for brevity in I, extend our work in several directions, 
and provide proofs for various conclusions drawn in I. Among these is a proof of a theorem 
(see Section 5) describing conditions under which certain random growth processes avoid 
intersection in d dimensions; as such, its utility extends beyond its particular use in this 
paper. 

The outline of the paper is as follows: In Section 2, we define our model. In Section 3 we 
provide a simple algorithm for finding the ground state configuration for any volume with a 
specified boundary condition, and discuss some of its properties. In Section 4 we demonstrate 
a second algorithm for finding the ground state, thereby mapping our problem onto invasion 
percolation and showing how the ground state multiplicity problem in our model is equivalent 
to the question of nonintersection of invasion regions in invasion percolation. In Section 5 we 
provide a heuristic argument for the dimension-dependence of the above questions, make part 
of that argument rigorous by relating it to the square-summability of a connectivity function 
in invasion percolation, and discuss the dimensionality of the invasion region as a function of 
space dimensionality. This yields the previously mentioned transition at eight dimensions. In 
Section 6 we discuss the role of frustration, and point out an important difference between the 
nature of multiple ground states in spin glasses and that in random ferromagnets. Finally, in 
Section 7 we summarize our results and briefly discuss their application to some dynamical 
problems of random walks on rugged landscapes [pi], [22 . 



2 Model 

We work on a cubic lattice in d dimensions, i.e., of sites x G Z d and edges connecting 
nearest-neighbor sites only. The Hamiltonian is of the standard Edwards-Anderson Ising 
form given in Eq. (fj). The difference with the usual EA Ising model is in the coupling 
distribution, which now depends on the system size. That is, we apply a nonlinear scaling 
to the couplings which "spreads them out" so much that each coupling magnitude exists on 
its own scale — more precisely, for large enough system size each coupling will have at least 
twice the magnitude of the next smaller coupling. 

We therefore consider a region A^, which is an L d cube centered at the origin. We 
achieve the required condition on the couplings by separating their signs and magnitudes 
in the following manner. Let e xy be a set of independent identically distributed (i.i.d.) 
symmetric ±1 valued random variables, and let K xy be a set of i.i.d. continuous random 
variables {e.g., uniform on [0, 1]). The e X yS and K xy y s are defined on a common probability 
space (Q, JF, P) and are independent of each other. A point u in Q may be thought of as a 
particular realization of all e^'s and K X y§. Then we set 

Jxy = J xy * = CL^xy G V j (2) 

where cl is a linear scaling factor which plays no role in ground state selection and where 
the nonlinear scaling factor A^ is chosen to diverge fast enough as L — > oo to ensure that 
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(with probability one) for all large L, each \Jjff\ in A^ is larger than at least twice the next 
smaller one. 

To see that such a choice of is possible, note that for any distinct pair of edges, the 
function, 

g(X) = P(~ < e - XK *«/e- XK *V < 2) (3) 

— P( \K X y — K X 'y'\ < IJ y), 

tends to zero as A —>■ oo because K xy and K x i y i are independent continuous random variables. 
The probability that our desired condition on the |J^|'s in is not satisfied is bounded 
by 

E E S(A (L) ) = 0(iV))), (4) 

<xy> < x 'yi> 

where the sums are over edges in A^. If we choose \^ so that g(X^) = 0(L~( 2a!+1+<: )) for 
some e > 0, then the sum of (4) over L is finite, so by the Borel-Cantelli lemma it follows 
that (with probability one) our desired condition will be valid for L > some finite L*(uj). 
For example, if the K X yS are uniform on [0, 1], then g{X) = O(j) and so A^ L ^ > L^ 2d+1+ ^ is 
a sufficiently fast divergence. 

It may be helpful to note that the couplings (2) can be expressed in a simpler form if 
we choose each A^ to be an odd integer. Letting J xy denote the L-independent coupling 
e xy exp(—K xy ), Eq. (2) then becomes 

4? =C L (j X yf L \ (5) 

Any continuous symmetric distribution for the J X yS (such as Gaussian) is possible. 

We will show in the next section that we have constructed a model whose ground state 
for large L can be found with a type of greedy algorithm. While this may seem to suggest 
that interesting ground state behavior cannot occur, we shall see that, surprisingly, this is 
not the case. We point out here, however, that the model has no interesting behavior at 
nonzero temperature. Its use is only as a means of studying ground state structure. 



3 Greedy Ground State Algorithm 

We begin with a formal definition of infinite-volume ground states for a specific coupling 
realization u. We first consider the finite volume A^ and for the moment assume some fixed, 
boundary condition, a, on 8Al, the boundary of A^. We define cr£ to be a ground state on 
Al (with boundary condition a) if it minimizes Hl, where Tii is the Hamiltonian of Eq. (H) 
with the summation confined to couplings in A^ (including those between A^ and 8Al). The 
set of all infinite volume ground states given u is then the set of all subsequence limits as 
L — » oo of o~2 with some &l (i.e., the set of all possible limits with all possible boundary 
conditions). 

For the volume A^ with a fixed boundary condition on d Al, we now describe an algorithm 
for finding the ground state a* L - (Non- fixed boundary conditions, such as free and periodic, 
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will be discussed in Section |6].) To do this we rank order the couplings in the following 
manner: the coupling with largest magnitude (corresponding to the smallest value of K xy ) 
in Al will be said to have rank one; the next larger coupling will have rank two, and so on. 
Now select the coupling with rank one; i.e., that with the smallest value of K xy , and choose 
the spins a x and o y on its endpoints such that a x e xy a y > 0, i.e., the coupling is satisfied. 
Then find the coupling of rank two, and choose the spins on its endpoints to satisfy it. Repeat 
this procedure, unless a resulting closed loop (or path connecting two boundary sites) with 
previously satisfied couplings forbids it. When that happens, simply proceed to the coupling 
next in order, and continue until every coupling has been tested. 

Note that until a cluster of spins (connected by tested edges) reaches the boundary, 
only the relative orientations of the spins in the cluster are known. With a fixed boundary 
condition, the sign of each spin in a cluster will be determined as soon as it connects to the 
boundary. 

It is not hard to see that this algorithm always provides the ground state for any in 
which every coupling magnitude is greater than the sum of all those of lower order in A^. 
This will be the case for all large L (with probability one) because each coupling magnitude 
is greater than at least twice that of the next smaller one. To see then that in the ground 
state a given coupling must be satisfied providing this doesn't violate the previously satisfied 
couplings of higher rank (or the boundary conditions), consider the clusters formed by the 
previously satisfied couplings. Under the proviso, the two endpoints of the given coupling 
must belong to distinct clusters, at least one of which does not touch OAl. If the given 
coupling were not satisfied in a spin configuration, then flipping all the spins in that cluster 
(the one not touching dA£) would lower the energy and so the spin configuration would not 
be a ground state. 

The algorithm outlined here is easily recognizable as simply a version of the greedy 
algorithm — in effect, we have invented a spin glass model whose exact ground states can be 
found via the greedy algorithm (which for most models is generally a relatively poor algorithm 
for finding ground state configurations or energies). It is not at all clear at this point, however, 
that the procedure that we outline, when repeated for ever-increasing volumes, will have a 
natural infinite-volume limit — but we will show in the next section that this is in fact the 
case. Before we do that, however, we explore some of the properties of our model in light of 
the ground state algorithm just described. 

3.1 Statement of the Problem 

The question of whether this model has multiple infinite volume ground state pairs is equiv- 
alent to whether, as L —>■ oo, a change in boundary conditions can change a fixed coupling 
deep in the interior from being satisfied to unsatisfied, or vice-versa. 

We therefore ask whether any bond < X1X2 > exists with the following property: for all 
large L, before any path of satisfied couplings joining x\ and 22 within A^ (i.e., not touching 
8Al) is formed according to the greedy algorithm described above, there already exist two 
disjoint paths, one joining x\ to the boundary, and the other joining x<i to the boundary. If 
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such a bond exists, then whether its coupling is satisfied or unsatisfied will be determined 
by the boundary conditions (for all large L). 

When the boundary of is sufficiently far from some fixed interior region R, it may be 
that no such coupling exists within R — each coupling in R is either itself tested before two 
such disjoint paths can be found, or else its endpoints are first connected via some path (not 
touching dA L ) of previously tested couplings. If that is the case for every finite R, then only 
a single pair of spin-flip-related ground states exists in the thermodynamic limit. Otherwise, 
the system possesses multiple pairs of ground states. 

3.2 "Always Satisfied" Bonds 

One can distinguish between two kinds of satisfied bonds in our model, for a given coupling 
realization u: there is the kind which is satisfied in uo but which would become unsatisfied 
in u/, which is simply uj with the sign (i.e., the e xy ) of that particular bond reversed. That 
is, whether this bond is satisfied depends on the specific sign of its corresponding e xy . The 
second set of satisfied bonds are those which are satisfied regardless of their sign; we will call 
them SI bonds (a precise definition is given below). It is easy to see, for example, that the 
coupling of highest rank is SI. In fact, each site in A^ has at least one SI bond attached 
to it, namely that whose rank is the highest of all bonds which connect to that site. We 
note that SI bonds appear in the ordinary EA model also (e.g., any bond whose magnitude 
is greater than the sum of the magnitudes of the adjoining bonds at either of its ends) but 
they don't appear to play the crucial role that they play in our model. 

Because SI bonds play an important role in what follows, we will devote some space 
to study their properties here. In particular, it is they which determine the ground state 
structure in our model. We begin with a precise definition of SI bonds for a given A^. 

Definition . A bond will be denoted SI if the following is true: its rank must be greater 
(i.e., its coupling must be of larger magnitude) than at least one coupling in any path (not 
using that bond) connecting its endpoints. In this definition, we treat all boundary points 
as automatically connected so that the union of a path from x to 8Al and a path from y to 
dAi is considered to connect x and y. We remark that this notion of connection within A^ 
is a consequence of our dealing with a fixed boundary condition a; when we treat boundary 
conditions such as free or periodic (in Section 6), the notion of connection will be modified 
accordingly. 

According to the definition, an SI coupling is chosen to be satisfied by the greedy algo- 
rithm before any other path of similarly chosen bonds connects its endpoints. It is apparent 
that it is the SI couplings which determine the ground state configuration. Satisfied bonds 
which are not SI (call them S2) play no role in determining any part of the ground state 
spin configuration. 

We now present some properties of SI bonds which will be useful later. 
1) The set of all SI bonds forms a union of trees. This claim is obvious from the definition 
of SI bonds. 
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2) The set of all SI bonds spans the set of all sites in A L , i.e., every site belongs to at least 
one SI tree; and furthermore every SI tree touches the boundary of A^. 

This second claim comes in two parts: the first has already been shown above by explicit 
construction. The second part is easily shown by contradiction: suppose a given SI tree 
"dies" before reaching the boundary of A^. Consider all edges which connect a point in this 
tree to a point not in it. The coupling of highest rank within this set must be SI. Therefore, 
all trees formed of SI couplings reach the boundary. 

It may happen, for a given u, d, and A^, that the SI couplings form either a single tree 
or a union of disjoint trees. Note that this tree partition of A^ is the same for all boundary 
conditions a. Within each tree the relative sign of the spins is fixed by the SI couplings; the 
overall sign for each tree is determined by a. We can now address the question posed earlier 
- different boundary conditions can give rise to different infinite volume ground states if and 
only if there exist fixed neighboring sites which belong to disjoint trees of A L for arbitrary 
large L's. 

We are left, however, with several important questions — how is the tree partition for a 
particular L related to that for some V > LI More specifically, how can one be sure that 
the procedure we have proposed has a natural infinite-volume limit? In order to answer 
these questions, we present an alternative algorithm in the next section, which will provide 
a mapping to invasion percolation. 

4 Invasion Percolation Algorithm 

Before describing our alternative algorithm for obtaining the ground states, let us note that 
although the absolute ranks of the K xy 's depend on L, the K X yS themselves, and hence their 
relative ordering by rank, does not change with L. This will allow us to analyze the L — > oo 
limit of our algorithm. 

We begin by defining (in all of Z d ) and for a given u (and hence a fixed relative ordering 
of all the K xy (ujys), a growing sequence of trees T (u), Ti(u), T 2 (u), . . ., starting from u e Z d , 
with T n {u) containing n edges and n + 1 sites (including u). T (u) consists of u alone, and 
in general T n+ i(u) is obtained from T n (u) by considering all edges from sites in T n {u) to new 
sites and adjoining the edge e n+ i(u) =< x n ,x n+ i > (and new site x n+ i) with the smallest 
value of K xy . This procedure is identical to that employed in invasion percolation on Z d 
(except that we include only edges which connect to new sites). 

For a given L and u G A L , let N L (u) denote the smallest n such that T n {u) touches dA L . 
The crucial point is that when every coupling magnitude in Al is greater than the sum of all 
those of lower order, then for any boundary condition a on OAl (and any choice of e XJ ,'s), in 
the ground state a* (a), every coupling in T^ L \u) = T Nl ^(u) must be satisfied. To see this, 
note that if J XnXn+1 were not satisfied (here n + 1 < N L {u)) in a given spin configuration, 
then flipping all the spins in T n {u) would lower the energy because J XnXn+1 is the coupling 
on the boundary of T n (u) of largest magnitude, er* is then determined by the tree T^ L \u), 
the coupling signs e xy on that tree and the boundary condition a x at the boundary site x 
touched by that tree. 
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It is clear from the preceding discussion (and that of Section 3) that every bond in 
T( L '(u) (for every u) is an SI bond. Since every T^ L \u) touches <9A^ and the union of all 
these edges (for all it's in A^) clearly touches every site in A^, it must be that this union 
is the same tree partition of A^ as obtained in Section 3 from the union of all 5*1 bonds. 
On the other hand, it is clear from the last paragraph, that (for u,v & A L ) the relative sign 
cr*(a)o"*(o") is the same for all choices of a if and only if T^ L \u) and T^ L \v) are nondisjoint. 
Furthermore, one has the following dichotomy concerning the infinite volume invasion trees, 
Too( u ) — ^ m n^ooT n (u) = lim^oo T^ L \u): If T^u) and T^iv) are nondisjoint, then T^iu) 
and T^ L \v) are nondisjoint for all large L; if T^u) and T^v) are disjoint, then T^ L \u) 
and T^ L \v) are disjoint for all L (such that u,v E Al). We are thus led to the following 
conclusions, concerning the trees T^u) and their union, 

Foo = U neZ dT 00 (-u), 

which we call the invasion forest (note that both the trees and the forest depend only on the 
K X y J s and not the e xy 's): 

(1) In every infinite volume ground state, every coupling in is satisfied. 

(2) is either a single (infinite) tree or else a union of M > 2 distinct (infinite) trees; in 
either case it spans all of Z d . 

(3) The former case happens if for every u, v, the trees T^iu) and T^iy) intersect. In this 

case there is a single infinite volume ground state pair. 

(4) The latter case happens if T^u) and T^v) are disjoint for some u,v. In this case the 

number of ground state pairs is 2 Ar_1 and so is uncountable if M is infinite. 

In the next section we will discuss the dependence of the value of jV, and hence of the 
ground state multiplicity, on the spatial dimension d. 

5 Nonintersection in invasion percolation 

We have mapped the problem of multiplicity of ground states in our model to that of whether 
invasion percolation has nonintersecting invasion regions. We can already answer the ques- 
tion of multiplicity of states of our model in two dimensions. Because it is known that for 
d = 2 invasion percolation the trees T^^u) and T^v) always intersect (in fact are the same 
modulo finitely many sites) |§, it follows that our model has only a single pair of ground 
states in two dimensions. Whether any two such trees in higher dimensions must intersect 
is an interesting problem in invasion percolation, which we now consider. 

To proceed (mostly nonrigorously) , we use a well-known feature of invasion percolation: 
that the invaded region asymptotically approaches the so-called incipient infinite cluster (i.e., 
at the critical percolation probability p c ) in the independent bond percolation problem on 
the same lattice |28[]. The fractal dimension D of the incipient cluster in the independent 
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bond problem on the d- dimensional cubic lattice is known from both numerical studies and 
scaling arguments; in particular, D is dimension-dependent (increasing with d) below six 



dimensions, but D = 4 for d > 6 28 



The following heuristic argument might then provide an intuitive picture of our model's 
behavior. Consider the infinite- volume invasion trees T^u) and T^v) introduced in the last 
section. If each has a fractal dimension less than d/2, the probability that they will "miss" 
each other is greater than zero; if it is greater than d/2, they will intersect with probability 
one. 

This suggests that if the fractal dimension Di of T^u) is equal to D, then the critical 
dimension of our model is eight. Below eight dimensions invasion regions should always 
intersect, and hence there would be only one pair of ground states in our spin glass model; 
above that there should be an infinite number. 

To make this line of reasoning a bit more precise, let us define a pair connectedness 
function G(y — x) as the probability that the site y G T^x). We can then define Di by the 
relation |L1], |8f 

G(y-x)~l/\\y-x\\ d - D > (6) 

as || y — x || — > oo. Note that summing Eq. @ over all y in a box of side length L centered at 
x yields L Di as the order of the (mean) number of sites in that box which belong to T^x). 

Our main task is then to determine Di as a function of space dimension d. First, however, 
we provide a precise statement of a condition for nonintersection of invasion trees, which also 
provides rigorous justification for part of the heuristic argument presented above. 

Theorem . If J2 x ez d G(x) 2 < oo, then (with probability one) there are infinitely many 
(random) sites X\, x%, . . . such that T^Xi) fl T^Xj) = 0. 

Remark . Although we state this theorem in the context of invasion percolation, the 
proof we will now present shows that it remains valid for a fairly general class of random 
growth processes in place of T n (x). The ingredients of the proof are (statistical) translation 
and reflection invariance and the fact that the events {T n (x) = A} and {T n (y) = B} are 
independent as long as A and B are separated in Z d by some fixed distance. This fact is in 
turn a consequence of the "local dependence" of events like {T n (x) = A} on the underlying 
K xy variables (and the mutual independence of those variables). 

Proof . Let A n denote the event that there exist some n sites, xi, . . . ,x n with T^Xj) fl 
T<x>{xj) = for 1 < i < j < n. If we can show that for every e > and every n, P(A n ) > 1 — e, 
then P(A n ) = 1 for each n and the desired result follows by letting n —> oo. The desired lower 
bound on P(A n ) would itself be a consequence of showing that PiT^x) fl Too(?/) = 0) — > 1 
as \\x — y\\ — > oo. To see this, pick deterministic y±, . . . ,y n with \\yt — yj\\ large enough, for 
i ± j, so that P(T 00 (y i ) n T^y,) ^ 0) < e/Q; then 

P{A n ) > PiTooiyi) n Toc(%-) = for 1 < % < j < n) 
> 1- E P(T oo (y i )nT oo ( % -)^0) 

l<i<j<n 
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It remains to obtain a suitable lower bound for P(T oc (a;) nT^y) = 0) when ||sc — y\\ is large. 

Denote by p (T 00 (x),T 00 (y)) the minimum Euclidean distance between some site in T^x) 
and some site in T^y). Furthermore, let T^(y) denote an invasion region constructed using 
a completely independent duplicate set of variables {K' xy } (so that T^x) and T'^y) are 
independent random trees). An elementary but crucial observation is that 

PiT^^T^y) = 0) > P(p(T 00 (x),T 00 (y)) > 1) = P {pp^x),^)) > 1) , (8) 

where the inequality is trivial and the equality follows from a fairly standard type of argu- 
ment (see, e.g., Ref. 0) given in the next paragraph which, roughly speaking, leads to the 
conclusion that T^x) and T^y) are independent, conditional on p (T ao (x),T ao (y)) > 1. 

To explain more concretely the equality in Eq. (|3|), let us note that for any possible 
configuration A of T n (x), the event that T n (x) = A depends only on the K X yS with < xy >G 
S(A), where S(A) denotes the set of nearest neighbor edges which touch either one or two 
vertices of A. If B is a possible configuration for T n (y) with p(A,B) > 1 (or equivalently 
with £{A) n£(B) = 0) then the events {T n (x) = A} and {T n (y) = B} are independent since 
they depend on disjoint sets of the independent K X yS. Hence 

P(p(T n (x),T n (y))>l) = £ P(T n (x)=A,T n (y) = B) 

A,B:p(A,B)>l 

E P(T n (x)=A)(T n (y) = B) 

A,B:p(A,B)>l 

£ P(T n (x)=AX(y) = B) 

A,B:p(A,B)>l 

= P(p(T n (x)X(y))>l) ■ (9) 



Letting n — > 00 gives the desired equality. 

Let I z denote the event that z G T^x) and p{z ) T' oa {y)) < 1. Then p(T 00 (x),T^ x> (y)) < 1 
if and only if I z occurs for some z and hence 

P(p(T 00 (x),T^ (y)) >l) = l-P(U zeZ J z ) 

zez d 

= i-E p (*e T »(^(p(«X(y))<i) 

zez d 

> l-Y,G{z-x) E G(z'-y). (10) 

zez d ||z'-z||<i 

By using the reflection invariance of G(x), the last expression can be rewritten as 

!- E ( E + J . (11) 

IM|<1 \zez d / 
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To complete the proof, it clearly suffices to show that 



(G * G) (x) 



zez d 



z)G(z) 



as x 



•x- 



But 



where 



(G*G)(x) 



G(k) 



—ik-x 



dk 



G(k) = (2n)- d/2 J2 G(x)e' 
xez d 

Furthermore, G(k) is real since G(x) = G{—x) and so 



k-x 



(12) 

(13) 
(14) 



G{k) 



dk 



G(k)\ dk= \G{x)\ 2 < oo ; 
xez d 



(15) 



thus 



(G*G) 



G(k) 



x 



G L 1 ([— 7r, 7i] d , dk) and so by Eq. ([13]) and the Riemann-Lebesgue lemma, 
► as llxll — > oo, as desired, o 



Using this result and Eq. (|6]), it follows that the condition Di < d/2 is sufficient for 
nonintersection of invasion regions in d dimensions, and correspondingly, for our model to 
have an uncountable number of ground states. We now consider this question. 

Monte Carlo simulations of invasion percolation on square and simple cubic lattices P§| 
provide strong evidence that Di = D in dimensions two and three (D = 91/48 and D = 2.53 
In higher dimensions, less is known. There does exist, however, an exact solution for 



invasion percolation on a Cayley tree [23 



One can deduce the fractal dimension Di of the invasion region using two different mea- 
sures. One of these is to compute the radius of gyration (i.e., the root of the mean square 
cluster radius) R after the invasion process has completed n steps. For the Cayley tree, it 



n 



1/4 



consistent with 



was found using Monte Carlo simultations |23| that, as n — ► oo, R 
A = 4. 

The second measure uses the exact solution mentioned above and is much closer to 
Eq. (|6|). The shape function, S^, which is the mean number of invaded sites on level m of 
the Cayley tree for an invasion of n steps, is computed. By analyzing Eq. (9) in Ref. P3| 
(valid for the simplest Cayley tree, i.e., with coordination number 3) in the limits /3 — > 1 and 
a — ► 1 (in that order), we find that S!^ is proportional to m (with logarithmic corrections) 
as m —>■ oo. The total number of sites invaded up to level m thus scales as m 2 . The usual 
measure of distance on a Cayley tree places level m at distance y/m from the origin, leading 
again to Di = 4. 

We therefore conclude (nonrigorously) that Di = 4 is an upper bound for the fractal 
dimension of an invasion tree on a lattice in finite dimension. Given that we expect Di = D 
in any dimension, and given the known values of D, we conclude that the critical dimension 
in our problem is eight. 
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6 The Role of Frustration 



We have so far argued that our model has a transition in ground state multiplicity at eight 
dimensions. However, the nature of the "spin-glassiness" of our model, and in particular 
the role played by frustration, has not been clarified. In fact, suppose that one were to 
construct a model of a random ferromagnet using the Hamiltonian ([]]) and couplings (@) 
but with all e xy = +1. For this system, the ground state in any finite volume with specified 
boundary conditions can be found using the same greedy algorithm described in Section [3[ 
It is clear that once again many (infinite volume) ground states (in d > 8) can be generated 
with appropriate choices of fixed boundary conditions. (We recall that this is also the case 
for uniform ferromagnets in any dimension or ordinary random ferromagnets in d > 5). 

So how is the ground state structure of the above random ferromagnet different from 
that of our spin glass? The answer can be found through analyzing the behavior of each in 
the presence of spin- symmetric, coupling-independent boundary conditions, such as free or 
periodic. 

In an earlier paper, [Tj| the authors argued that multiplicity of ground states in the EA 
Ising spin glass should be associated with nonexistence of a single limiting Gibbs distribution, 
in the thermodynamic limit, for any coupling- independent boundary conditions. We find 
that the same association holds here. While the conclusion is fairly clear for any sequence 
of fixed boundary conditions, in which each boundary spin is assigned a definite value, the 
mechanism for spin-symmetric (e.g., free or periodic) boundary conditions is more subtle. 
Its investigation provides deeper insight into the nature of our spin glass model, and in 
particular the role played by frustration. 

We first discuss fixed boundary conditions. Let us now denote the always satisfied bonds 
in Al, as defined in Sec. [|, as SV" bonds, where the w superscript (for "wired") is to 
remind us of the convention that all boundary points are treated as automatically connected. 
(This is analogous to wired boundary conditions in the Fortuin-Kasteleyn random cluster 
representation of Ising and Potts models Jl|.) Let us denote by F™ the union of all the SV 
bonds, or equivalently the union over u E Al of the invasion trees T^ L \u) (stopped when 
they touch d Al)- F^ depends on the K X yS while the overall "sign" of any individual tree in 
F™ (i.e., the sign of a single spin in that tree) is determined finally by the boundary spin (at 
the particular site where the tree touches 8Al) and the coupling sign e xy of the SV bond 
touching that boundary site. For d > 8, as L increases with a given sequence of fixed, 
coupling-independent boundary conditions, the sign of each tree will randomly flip, and no 
single limiting ground state is obtained in the thermodynamic limit. 

Indeed the size dependence is such that all of the uncountably many ground states will 
appear as limits along coupling-dependent subsequences of volumes. To see that, we first 
condition on all the K xy 's and let Ti,T2, ••• be a list (in some order) of all the distinct trees 
from Fqo. Let r]^ denote the sign of (some fixed a Xi in) T i} as determined by the e^'s and 
the boundary condition a^ L \ From the previous discussion it should be clear that (i) for 
each L, r]^ = : ^ 6 consists of independent random variables, equally likely to 

be +1 or —1, and (ii) for varying L, the t/^'s are independent. It follows (for almost every 
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coupling realization) that for every m and for every assignment 77^ of signs rjf 1 = ±1 for 
i G {l,...,m}, the collection {L : r]^ = fjf 1 ^ for % = l,...,m} is an infinite subsequence of 
L's. By a standard subsequence diagonalization argument, it follows that for every one of the 
uncountably many assignments Tj of signs T] i = ±1 for all % > 1, there exists a subsequence 
{Lj(rj) : j > 1} so that for all 2, 77} — > 77^ as j — > 00. 

Consider now the case of the finite cube C Z d , with c/ > 8, and with, say, free 
boundary conditions. (We briefly discuss periodic and antiperiodic b.c.'s later.) At first 
it might seem that the connection between ground states and invasion percolation is no 
longer valid, because for every finite volume there is only a single cluster of always satisfied 
couplings. This is because, unlike when the boundary spins are fixed, the procedure of 
satisfying couplings of successively smaller magnitude will be continued until all sites are 
connected. To explain this more fully, we define a bond < xy > with both sites x and y in 
Al to be an SI? bond exactly as in the definition of an SI bond in Sec. ^ except that a 
path is said to connect the points x and y only if the path stays entirely within A^, never 
touching the boundary. In this case, it is not hard to see that the set F[ of all Si/ bonds 
forms a single tree spanning all of A^. 

To clarify the relation between free b.c. ground states and the invasion forest F^, we 
begin with a simple but important observation, that every Sl w bond with both endpoints 
in A L (i.e., leaving out any Sl w bond touching dA L ) is also an Si/ bond. Let us denote 
by E L the set of all Si/ bonds which are not also Sl w bonds. These are the bonds which 
connect together the distinct trees of F™ to create the single spanning tree F[. Note that 
every such bond < xy > must have its two endpoints in distinct trees of F™. We now make 
a crucial claim: as L — > 00, the edges of E L will move out to infinity. More precisely stated, 
the intersection of with any fixed set of bonds (say, all bonds in Al with L fixed) will, 
for all large L, be the empty set. To see this, first observe that any such edge < xy > not 
moving out to infinity would have two properties: (a) There would be no path in Z d between 
x and y with all K values along the path strictly below K xy ; (b) T OQ (x) fl T^y) = 0. It was 
proved by Alexander [[| that such edges do not exist in any dimension. 

We will next argue that the relative sign between trees flips randomly, and so for our spin 
glass model there is no single limiting pair of ground states for a sequence of volumes (chosen 
independently of the couplings) with free boundary conditions. In fact, we will see that in 
our model, all of the uncountably many ground state pairs arise via coupling- dependent 
subsequences. This chaotic size dependence, i.e., absence of a limiting ground state for 
spin-symmetric, coupling-independent boundary conditions, is similar to that which would 
be found in the EA model were it to possess many ground state pairs [I9|. We argued in 



Ref. [nj that this is, in fact, the signature of multiple ground state pairs in spin glasses, 
and sets them apart from other systems which may also possess multiple ground state pairs, 
such as random ferromagnets. 

The argument for random flipping of relative signs between trees in the free b.c. case is 
like the one used above for flipping of absolute signs of trees in the fixed b.c. case, with a 
few modifications. Again, we condition on all K xy 's and let Ti,T 2 , ... be the distinct trees of 
Fqo. Now let 7)ij be the relative sign a Xi a Xj in the free b.c. ground state. This is determined 
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(finally) by the signs e xy of the bonds < xy > in E-. Since these bonds move off to infinity, 
we can pick a (sufficiently rapidly growing) subsequence Lk — ► oo so that the sets are 

disjoint from each other for different fc's. It follows that the relative signs rj[f between 
the first and i'ih trees will be independent as both i and k vary. Consequently, the same 
diagonalization procedure as in the fixed b.c. case can be applied to obtain the requisite 
further subsequences of L k . 

Finally, we note that for periodic (or antiperiodic) b.c.'s, Sl p bonds are defined by re- 
garding every boundary point as automatically connected to its periodic image point, but 
with no other automatic connections. Thus, like in the free b.c. case, every Sl w bond is also 
an Sl p bond, while the other Sl p bonds must move to infinity as L — > oo. This leads to the 
same conclusions as in the free b.c. case. 

Let us now contrast this picture with that of the random ferromagnetic version of our 
model. There, a sequence of coupling-independent spin-symmetric boundary conditions 
would yield a single pair of ground states as L — > oo — namely, all spins up and all spins 
down. This will be the case in any dimension. As with the usual models of ferromagnets and 
spin glasses, either of which might have many ground states, the difference from spin glasses 
is revealed most sharply through the use of coupling-independent, spin-symmetric boundary 
conditions. 

7 Conclusions 

In this paper we have expanded our original discussion in I, supplied additional proofs and 
arguments, and looked at our model in greater depth. In this section we briefly discuss 
some consequences of our studies, particularly what (if any) conclusions can be drawn about 
more realistic spin glass models, and how our approach can be used to study other topics, 
particularly certain dynamical problems. 

We first recap our main conclusions, based on both rigorous and nonrigorous arguments: 

• In invasion percolation on Z d , the invasion forest F^, has only one tree in dimensions 
less than eight, and infinitely many above eight. 

• In the spin-glass model introduced in I and discussed here, there exist only a single 
pair of ground states below eight dimensions, and infinitely many above eight. 

• Elucidation of the ground state structure exactly at eight dimensions cannot be ad- 
dressed by the techniques provided in this paper, and requires further study. 

• The random ferromagnetic verison of our model (all e xy = 1 in Eq. (Q)) also has 
a transition in ground state multiplicity at eight dimensions. However, exactly as in the 
Edwards-Anderson model, the special feature of spin-glass ground state multiplicity man- 
ifests itself in chaotic size dependence for spin- symmetric, coupling-independent boundary 
conditions, such as periodic or free. 

• The special properties of this model render it unsuitable for drawing firm conclusions 
about more realistic spin glass models. However, it suggests that a third conjecture (be- 
sides those of the Parisi picture and the scaling picture) should be entertained — that of 
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a transition in ground state multiplicity at some finite (and, obviously, greater than one) 
dimension. 

• Our study can resolve one perennially cloudy issue in the literature; namely, that of 
whether the joint presence of quenched disorder and frustration has any universal implica- 
tions for ground state structure or multiplicity The results presented here make it clear that 
no such a priori implications can be drawn. 

We close with a brief discussion of applications of our results to other problems. The 
form of the exponential expression in Eq. (Q) suggests that our techniques might apply in 
some way to the problem of a random walk in a strongly inhomogeneous environment. If 
K xy is thought of as an energy barrier between sites x and y, and A now represents inverse 
temperature, then exp(— \K xy ) may be interpreted as an Arrhenius factor proportional to 
the rate of making a transition from x to y at fixed temperature A -1 . Random K X yS, like 
those used throughout this paper, correspond to a rugged energy landscape. We can then 
prove that, in the limit as temperature gets small, the order in which sites are visited for the 
first time corresponds exactly to the invasion order in invasion percolation [pi]] . We further 
develop conclusions drawn from the theorems proved in Ref. |pT|| , and study the nature of 



broken ergodicity on a system whose state space resembles a rugged landscape, i.e, has many 



metastable states with random fixed barriers separating them |2Sq| . This leads to a surprising 
degree of emergent structure from what appears initially as a rather featureless landscape. 
The multiplicity of trees in the invasion forest for high dimensions plays an important role 



in that work. We are presently studying how our analysis here and in Refs. pi , |22j can be 
applied in other contexts. 
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